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ARTICLE INFO ABSTRACT 
Keywords: The prediction of ductile failure is crucial for the progress of metal forming industries. Therefore, 
Ductile fracture uncoupled fracture models have been proposed and continuously extended, trying to improve 


Uncoupled model 
Calibration procedure 
Number of tests 

Type of tests 


their forecasting abilities. However, without an assessment analysis of their predictive ability, it is 
difficult to establish the most suitable model, for describing the onset of failure under a wide 
range of loading conditions. Actually, a fair comparison of models’ predictive ability requires the 
use of an extensive experimental database and a similar calibration procedure. In this context, 
seven uncoupled fracture models are chosen and revisited to study the effect of the type and 
number of tests used for calibration on their predictive ability. Experimental datasets available in 
literature were considered for 2024-T351 aluminium alloy and for DP600 steel. For each material, 
the experimental dataset was split in two subsets. The first one includes tests that were carried out 
under plane stress conditions that were used in the calibration procedure. For the aluminium 
alloy, different groups of 6 tests, covering distinct stress states, were elicited from the experi- 
mental subset and applied in the calibration of the seven models. The remaining tests were used 
to assess the predictive ability and their sensibility regarding the group used in models’ cali- 
bration. Different trends of the fracture loci are observed, depending on the selected models and 
the tests used in their calibration. However, the analysis of the results indicates that a set of at 
least 6 experimental tests is appropriate to calibrate fracture models, providing that they enable 
covering a wide range for the stress triaxiality and Lode parameter. Finally, the conclusions 
provide some recommendations for the improvement of the calibration procedure. 


1. Introduction 


Thanks to the advent of digital computer capacities, numerical simulation codes have become fundamental tools in the develop- 
ment of metal forming processes, due to their ability to predict the material behaviour and occurrence of defects. However, the 
modelling of the ductile fracture is still a major challenge for the failure prediction during metal forming processes. Various approaches 
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Nomenclature 

eP Equivalent plastic strain 

€¢ Equivalent strain to fracture 

Eet Equivalent strain to fracture under uniaxial tension 

Ets Equivalent strain to fracture under pure shear 

a Equivalent strain to fracture under normalized uniaxial compression 
ZA Equivalent strain to fracture under normalized plane strain 

ah) Equivalent strain to fracture under normalized uniaxial tension 
Efpreaictea,  PYedicted fracture strain for the test i 


Experimental fracture strain for the test i 


Efexperimental, i 


0 Lode angle 

6 Normalized Lode angle 

é Normalized third invariant of the deviatoric stress tensor 
o Cauchy stress tensor 

Om Mean stress or hydrostatic stress 

On Normal stress 

01,02,03 Principal components of the Cauchy stress tensor 
o von Mises equivalent stress 

Of Equivalent stress to fracture 

Gur Hosford equivalent stress 

T Shear stress 

Tmax Maximum shear stress 

Nn Stress triaxiality 

(Wy, ----;W,) Material parameters 

A,n Hardening law parameters 

a,b,c Fracture parameters’ for the H-C model 


C1,C2,C3 Fracture parameters’ for the MMC model 

C1,C2,c3 Fracture parameters’ for the L-H model 

C1,C2,C3,C4,m Fracture parameters’ for the X-W model 

C;,C2,C3 Fracture parameters’ for the Q-K model 

Cc Height of cut-off value of stress triaxiality for L-H model 

c3*,c5,cy MMC fracture parameters’ that controls the dependence of the Lode angle parameter and the asymmetry of fracture 
locus, respectively 

D Damage indicator 

D,,D2,D3,D4 Fracture parameters’ for the M-W model 

D,,D2,D3,D4,Ds,D¢ Fracture parameters’ for the B-W model 

Err Relative error function 

fi.fo.fz Lode angle parameter functions presented in H-C model 

f(n,L) | Dependency function on Lode parameter for the L-H model 

&i[n,L] | Dependency function on stress triaxiality and Lode parameter for Q-K, 2020 model 


gL] Dependency function on Lode parameter for Q-K model 
i Number of test 

J2,J3 Second and third invariant of the deviatoric stress tensor 
E Lode parameter 

N Number of tests used for calibration 

$1,52,53  Deviatoric principal stresses 

s Deviatoric stress tensor 

() Macaulay brackets 


have been adopted to improve the understanding of the mechanisms involved in ductile fracture, by combining experimental and 
mechanical testing ([1,2]); microscopic observations ([3,4,5]) and numerical analysis ({6,7]). Nowadays, it is consensual that the 
nucleation, growth and coalescence of micro-voids play the major role at the microscopic level, influencing shear band formation 
({8,9,10,11]). Therefore, several damage models have been proposed based on these ductile fracture mechanisms, which are 
commonly classified into coupled and uncoupled models. 

The first category of ductile fracture models comprises coupled models, where the damage accumulation is inherent to the 
constitutive model, i.e., the stress tensor is coupled with the damage variable that affects the dissipation potential of the plastic 
deformation. The most widely known coupled ductile fracture model is the one commonly referred as Gurson-Tvergaard and 
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Needleman (GTN) [12]. It is based on the Gurson’s model [13], which assumes the damaged material as a porous media (porous 
plasticity theory) and, consequently, introduces the concept of void volume fraction. However, the void growth assumption depending 
only on the stress triaxiality is presumably insufficient to judiciously describe the ductile fracture behaviour. Thus, Tvergaard and 
Needleman extended this model to account for nucleation and coalescence mechanisms and included also some parameters to improve 
the correlation with results obtained from unit cell computations [12]. Although, Gurson-like models can accurately describe the 
damage evolution for several applications they also revealed some drawbacks, namely: (i) no void shearing effect is included [14]; and 
(ii) due to the large number of parameters, its calibration can be a costly and laborious task. Therefore, to overcome these deficiencies, 
several extensions have been proposed, particularly to enable the successful prediction of failure under shear dominated stress states, 
by taking into account the sensitivity of the void volume fraction related to the Lode parameter ([15,4,16]). 

An alternative type of coupled models was developed based on the phenomenological continuum theory, commonly referred as 
continuum damage mechanics (CDM). In this framework, the dissipation potential was built based on thermodynamic principles and 
the equivalent stress is coupled with a damage variable to describe the degradation of the material. The most widely known CDM 
model is the Lemaitre [17]. Several extensions of the original Lemaitre’s model have also been proposed. However, these type of 
models often suffer from inaccurate predictions of strain localisation for pure shear stresses, (i.e., limited range of stress triaxiality) and 
particularly, under complex strain paths [18]. Some modifications have been introduced to overcome these limitations (see e.g., [19]). 
Nonetheless, it should be mentioned that when using coupled models, the constitutive models that describe the plastic and fracture 
behaviour need to be calibrated simultaneously. In this context, several calibration strategies have been proposed based on destructive 
and non-destructive tests. However, it is difficult to determine a unique set of parameters, due to the interaction between the material 
plastic and fracture behaviours. Moreover, the evaluation of the micromechanical based variables is an onerous task, since the analysis 
of voids evolution requires the use of techniques such as, tomography, metallographic observations, acoustic emission or induction 
spectroscopy [20]. 

The second category of ductile fracture models encompasses the uncoupled models. These models have been increasingly devel- 
oped, since they present some advantages when compared with coupled models. Indeed, for uncoupled models, the damage accu- 
mulation is independent of the plastic dissipation potential. Moreover, their simple analytical forms, makes its calibration and 
implementation in numerical codes quite easy tasks ({21,22,10]). Besides, some works reported their good performance in failure 
prediction under complex loading paths ([18,23,24]). 

The uncoupled ductile fracture models have been developed based on experimental investigations at both macroscopic and 
microscopic scales. On one hand, from the macroscopic point of view, fracture results from the progressive degradation of material, 
since it gradually losses its stiffness and consequently, its strength. On the other hand, from the microscopic point of view, the material 
fracture is the consequence of: (i) nucleation, growth and coalescence of micro voids, under large stress triaxiality; (ii) the presence of 
shear bands, related to low stress triaxiality; and (iii) combined shear and void growth mechanism, attributed to medium range of 
stress triaxiality. 

Cockcroft and Lathman, (C-L) model is among the first proposed uncoupled ductile fracture models [25]. This ductile fracture 
criterion considers that the main factor leading to material failure is the maximum principal tensile stress. It has been widely used for 
the prediction of fracture for bulk forming processes, such as forging [26] and in high-velocity perforation and penetration problems 
[27]. Nevertheless, its predictive fracture performance was shown only for the region of low and negative stress triaxiality 
({28,29,30]); thus, it is restricted to this range [26]. This limitation is certainly related to its low flexibility (it involves only one damage 
parameter) ({26,28,31]). 

Wierzbicki and co-workers assumed that when damage occurs in metallic materials, the accumulated equivalent plastic strain, 
expressed as a function of the stress state, reaches a critical value known as fracture strain [32]. It is known that the stress triaxiality 
influences the microscopic damage mechanisms, as previously mentioned. Thus, Mae and Wierzbicki proposed an uncoupled model 
(M-W) which only considers the dependence of fracture strain on the stress triaxiality and they evaluated its predictive ability using a 
combined experimental-numerical approach ([33,34]). Later, Wierzbicki and co-workers have shown that the fracture strain is limited 
within two bounding lines, corresponding to generalized shear or plane strain (lower bound) and to axisymmetric stress state (the 
normalized uniaxial tension or compression) (upper bound), respectively ([2,26]). They proposed the X-W model which defines the 
limits of the symmetric 3D fracture locus by a set of exponential functions. The two limits are defined by the effect of the normalized 
third invariant of the deviatoric stress tensor on fracture strain. However, the symmetric form of the fracture locus can limit the 
description of the material’s ductility [35]. Therefore, a 3D uncoupled fracture locus was proposed by Bai and Wierzbicki (B-W model), 
based on three bounding limits [32]. The lower bound is the same (generalized shear or plane strain), but the symmetry between the 
stress states corresponding to axial symmetry in deviatoric tension and compression is no longer imposed, leading to two distinct upper 
bounds. The B-W model takes into account the effect of hydrostatic and deviatoric stresses on the ductile fracture, based on the stress 
triaxiality and Lode angle parameter. Such kind of models were successfully used to predict ductile fracture ([32,36,37]), for various 
applications in automotive, aerospace and marine industries ([35,38,21]). Nevertheless, the main challenge for these models was their 
ability to predict failure under non-proportional and complex loading conditions, for different classes of metallic materials. 

Another research trend was followed for developing uncoupled ductile fracture models that is based on modified versions of 
fracture models. One of the most deep-rooted fracture model is the Mohr-Coulomb (M-C) [39], which has been extensively applied to 
predict fracture of non-metallic materials, such as soil and rocks (brittle materials, when compared to metals and their alloys). Since 
the mechanisms involved in the fracture behaviour of brittle materials are very different from those of ductile materials, the appli- 
cation of the M-C model was extended to metallic materials by introducing the influence of both the hydrostatic pressure and the Lode 
angle parameter on the fracture behaviour [40]. Bai and Wierzbicki successfully modified the M-C stress-based fracture criterion into a 
3D fracture locus (MMC model) [40]. The performance assessment of the MMC model was achieved by its application for a diversity of 
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metallic materials, as reported in ([41,42,43,31,44,40,9,26,45,46]). 

Despite the great success of the MMC model, it has been shown that it is consistently more suitable for failure prediction in shear- 
dominated stress states, rather than in ones controlled by void-growth and coalescence mechanisms [40]. Therefore, further modi- 
fications were brought to this model to improve its prediction accuracy for a wider range of stress states ([47,48]). Essentially, Tresca 
equivalent stress was substituted by the Hosford one, leading to the so-called Hosford-Coulomb (H-C) model [47]. The H-C model 
predicts fracture when the linear combination of the normal stress and the Hosford equivalent stress reaches a critical value. According 
to Ha et al. [49], the H-C model exhibits an improved performance for the low positive range of stress triaxiality in comparison with 
MMC model. It was applied to predict the ductile fracture initiation for several metallic materials, as reported in 
([50,51,24,47,48,52,23]). 

Additional uncoupled models were developed based on microscopic damage accumulation mechanisms, such as the nucleation, 
growth and shear of voids. Lou and Huh (L-H) model is as an example [29], which is characterized by the incorporation of a cut-off 
value, for which no fracture occurs, when the stress triaxiality is less than (-1/3). This value is consistent with the observations reported 
in [53]. The L-H model was implemented in Abaqus/Explicit solver to predict the fracture behaviour of DP780 steel [29]. The results 
showed that for non-proportional loading conditions, the cut-off value could occur for different stress triaxiality values. Thus, this 
model was further modified by introducing a variable which takes into account the material microstructure effect on the stress 
triaxiality cut-off value [54]. In addition, the influence of the anisotropic behaviour of the material on ductile fracture was also 
included in the model [55]. L-H model was successfully applied to predict fracture for isotropic ({29,10]) and anisotropic materials 
({56,55,57]). The results showed good fracture predictions under different loading conditions, such as tension and pure shear stress 
states [56]. 

Recently, Quach et al. proposed an uncoupled model to describe the anisotropic fracture behaviour of metallic sheets ({9,58]). The 
results showed that the Q-K model enables the prediction of fracture in a square cup drawing test performed with an 6014-T4 
aluminium alloy [9]. 

The above-mentioned uncoupled models were used to solve engineering problems involving ductile fracture prediction, including 
the description of localization phenomena observed under complex loading conditions or/and for non-proportional loading paths. 
Consequently, it is evident that significant research efforts have been deployed to reveal the inherent mechanisms responsible for the 
ductile fracture behaviour of metals and their modelling. Despite the large number of available ductile fracture models, the selection of 
the appropriate one to accurately describe the fracture behaviour of metals is still a challenging task. The models’ calibration results 
depend on the material, the number and the type of tests, which are related to the stress state range. Therefore, it is of paramount 
importance to study its effects, using several combinations of experimental tests, on the models’ ability to reproduce experimental 
observations. 

In this work, a calibration procedure is used to assess the effect of the set of experimental data selected on the parameter identi- 
fication and on the fracture prediction ability of seven uncoupled ductile fracture models, as such thorough analysis has not been 
reported elsewhere. In this context, the same calibration procedure is applied to several sets of experimental tests, defining different 
identification strategies for the 2024-T351 aluminium alloy. One of the calibration strategies is then applied to the DP600 steel to 
evaluate to what extend the drawn conclusions are valid for other materials. 

This manuscript is organized as follows. The next section presents the variables used to describe the stress state in uncoupled 
fracture models. A brief description of the seven uncoupled fracture models is presented in section 3. The available sets of experimental 
results for both materials are depicted in section 4, as well as the procedure used to identify the models’ ductile fracture parameters. A 
detailed discussion of the results is presented in section 5 and the main conclusions are drawn in the last section- 


2. Description of the stress state 


Ductile fracture mechanisms are related to the loading conditions that the material is undergoing, i.e., to the stress state, which can 
be defined by the Cauchy stress tensor, o. Regarding the ductile fracture mechanics, the stress state is usually described by two 
dimensionless parameters: the stress triaxiality 7 and the normalized Lode angle 6. These parameters can be defined based on the 
second and third invariants of the deviatoric stress tensor, s, expressed as follows: 


Jn = ss : 8, J; = det(s) (1) 


$= 6-01 ==(0, + 0. +03) (2) 


Wl 


where 0}, 62 and o3 are the principal stresses of o, such that 0 >02203. The von Mises equivalent stress can be expressed as follows: 


= 55 :s= V3h = v3 [(o1 02) + (0 — 03) + (65 —o1)'| (3) 


The stress triaxiality is defined as the ratio between the hydrostatic stress and the von Mises equivalent stress: 


(4) 
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The stress triaxiality is considered as a fundamental factor of ductile fracture for metallic materials. In fact, the hydrostatic stress is the 
driving force of the nucleation and growth of voids. On one hand, high hydrostatic stress intensity leads to an increase in the voids’ 
sizes. On the other hand, the deviatoric stresses are responsible for the change of the voids’ geometry. In this context, it is also 
important to define the normalized third invariant of the deviatoric stress tensor, as follows: 


21, 3V345 
~ 2 aR? 


(5) 


which characterizes the position of the second principal stress o2, with respect to the maximum o, and minimum 03 principal stresses, 
such that it is bounded between —1 and 1. This normalized parameter ¢ is related to the Lode angle 0 through the following relation: 


& = cos(30) (6) 


such that the range for the Lode angle 0 is 0<0<3. Besides, it is also common to adopt the normalized Lode angle 6, defined by: 


= 60 2 = 
é=1 == 1 ~ parcos(¢), where —1<6<1 (7) 


which is also called Lode angle parameter. 6 = 0 corresponds to generalized shear or plane strain, while 96 = —1 and +1 correspond to 


axisymmetric compression and tension, respectively. Finally, it should be noted that the standard expression of the Lode parameter, 
denoted by L, is defined as follows: 


20) —0,— 03;  3tand— V3 


Ls ; 
0| — 03 tand+ /3” 


such that — 1<L<1 (8) 


In the particular case of plane stress condition (PSC), the hydrostatic pressure changes along the loading history, but the normalized 
Lode angle remains constant. Thus, the stress triaxiality and the normalized Lode angle parameter are related by the following 
expression [36]: 


= 2 27 (, 1 2. 2 
ae —ten<2 
d=1 7 2005( 5 a(n :) : for 3S; (9) 
Alternatively, the normalized third invariant of the deviatoric stress tensor is expressed by: 
1 = 27 ge 
‘ 1-0 f 1 
g cos (5 ( )) 5 a(n :) (10) 


The assumption of plane stress state is widely adopted in the calibration strategies of uncoupled fracture models, because it 
simplifies the procedure due to the representation of the fracture locus in a 2D space (see Eq. (9)). Indeed, different types of engineering 
applications have been reported in literature, such as stamping of an anticlastic structure [60]), cross-die forming and cup drawing 
applications ([61,62]), incremental sheet metal forming [63], three-point bending ona hat assembly [52], laser-welded stiffened panel 
penetration tests [64], impact tests on flat and perforated plate [50], burring sheet metal [65], for which the PSC was assumed in the 
calibration procedure. Therefore, as far as we know, this assumption does not necessarily influence the accuracy of the predicted 
results. Moreover, the PSC is commonly accepted for thin walled metallic parts, including sheet metal forming processes. 


3. Overview of the selected ductile fracture models 


The seven uncoupled ductile fracture models selected represent a set of the ones that were proposed in, approximately, the last two 
decades. They all share a relatively low number of fracture parameters since they are phenomenological. In that sense, one could infer 
that the ones proposed in the first decade (i.e., the family of models developed by Wierzbicki et al. (X-W, 2005; B-W, 2007 and M-W, 
2008)) are more flexible, due to the slightly higher number of damage parameters than the one adopted in the more recent models (i.e., 
Q-K, 2020; H-C 2014; L-H, 2013, MMC, 2010). Nevertheless, different approaches were adopted, in particular, to integrate the material 
plastic behaviour in the description of the fracture surface. In this context, it is important to compare the predictive ability of more 
recent with older models, which were more widely investigated (e.g., X-W, 2007; B-W, 2007 and M-W, 2008). In this section, the 
selected models are briefly described. 


3.1. Xue and Wierzbicki, 2005 (X-W) model 


Bao and Wierzbicki [1] performed a series of experiments on a 2024-T351 aluminium alloy and confirmed that ductile fracture is a 
non-monotonic function of the stress triaxiality, as early suggested in [3] and [66]. However, the stress triaxiality is not sufficient to 
capture the influence of the multiaxial stress state on the material ductile failure behaviour. The important effect of the Lode angle 
parameter on damage prediction was shown in several studies (e.g., [67,43,26,40,5]). Both stress triaxiality and normalized third 
invariant of the deviatoric stress tensor were taken into account in the X-W model, which defines the fracture strain,é;, as follows: 


(1, é) Ceo" (Cyeo" C3e" 7) (1 — é m ye (1 1) 
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where C,e~°2"” and C3e~°*” are exponential functions that should be determined based on the equivalent strain to fracture under 


) (+ 


tek ee ee ee .. —(0 
uniaxial compression, a , uniaxial tension, é; a and plane strain, a ) such that: 


2° = a(n, € or 8 =0) = Ce" 
ay =8) =8(n,é or 0= +1) = Ce" 


The exponent m must be an even integer in order to enable the description of the drop of material ductility due to é by an elliptic 
function, with a symmetry plane in € = 0. Some authors suggested to select the even integer closest to 1/n, where n is the strain 
hardening exponent [2]. Therefore, the model has four damage parameters:C1, C2, C3 and C4. 


3.2. Bai and Wierzbicki, 2007 (B-W) model 


B-W model corresponds to an asymmetric ductile fracture model that neglects the effect of the strain hardening exponent on the 
fracture strain, unlike the previous model ([32,35]). The dependence of the fracture strain on 7 is described by an exponential function, 
whereas for é it is assumed as parabolic. This model was developed based on three bounding limits, i.e., three values of £, corre- 
sponding to axial symmetry in deviatoric tension, £ = + 1, plane strain or generalized shear, € = 0, and axial symmetry in deviatoric 


compression, € = —1, respectively. The fracture strain locus is defined by 
= 1 | = 
&(n,0) = 1 (Dye~P2" + Dse~P*") — Dye?" | + 5(Die™ — Dse~?*")6 + D3e?#" (12) 


In this case, the exponential functions should be determined based on the equivalent strain to fracture under uniaxial compression, 
( (+ 


a tension, é; , and plane strain, 2, such that: 
a) = &(n§ or 8 = —1) = Dye 

a =%(n, € or 0 = 0) = Dye?" 

ay = &(n,€ or 9= 1) = Die?" 


This model has six parameters: D,, D2, D3, D4, D; and Dg. The authors recommend the use of six types of specimen to carry out the 
experimental tests: tensile tests on smooth and notched round bars (0 = 1), plane strain test (9 = 0), flat grooved tension (0 = 0), 
upsetting test and cylinder upsetting test (9 = —1) [35]. 


3.3. Mae and Wierzbicki, 2008 (M-W) model 
This simple empirical ductile fracture model considers only the influence of the stress triaxiality, which is known to control the 


fracture mechanism of coalescence, growth and nucleation of voids, or/and shear mode ([33,34]). This model is expressed in three 
branches, as follows: 


ee! <n<0 
' T+3n' ae 
1 
B= Fret (Ex —Frs)(3N— 1), O<ns3 Of) 
_ e 1 
€¢ = Doe 37 + Day, 381 


where €;; and é;, are the effective fracture strain under uniaxial tension and shear, respectively. The four model’s damage parameters 
are D;, D2, D3 and D4. However, as can be seen from Eq. (13), Di = és and @ = Dze?s/3 + D4, which means that the parameters are 
not independent ([68,69]). M-W model was successfully applied in the prediction of ductile fracture for materials that experienced 
proportional loadings, for which the stress triaxiality remains constant [34]. 


3.4. Modified Mohr-Coulomb, 2010 (MMC) model 


Originally, the Mohr and Coulomb criterion [39] was developed for brittle materials, such as soil and rocks. It incorporates a linear 
relationship between the shear, zt and the normal oy, stress components and corresponds to the following relation: 


max(t + C10n)¢ = C2 (14) 
where, c and cz correspond to the model’s parameters, such that fracture occurs when the linear combination of the shear and normal 


stresses reaches the critical value cz. The main drawbacks of this model are the following: it can be only used under small plastic strains 
and it is pressure or stress triaxiality independent. Therefore, a modified version was proposed by Bai and Wierzbicki [40]. They 


M. Nouira et al. Engineering Fracture Mechanics 279 (2023) 108952 


introduced a hardening rule, which incorporates the dependence on the Lode angle parameter and stress triaxiality. The Mohr and 
Coulomb model is defined in a Cartesian coordinate system, based on the principal stresses (0 ,02,03). Hence, the first step consists on 
its transformation to the spherical coordinate system, defined by the equivalent stress, the stress triaxiality and the Lode angle (6,7,0). 
Assuming that the material is isotropic, the von Mises yield criterion can be adopted to determine the equivalent stress, as follows: 


-1 
yt Scos A)+e (n+ jana) (15) 


Further details on this transformation can be found in [40] and [70]. Considering that the material’s hardening behaviour can be 
described by a simple isotropic power-law. 


= A(e")" (16) 


Of = C2 


The equivalent strain at fracture can be expressed as: 
1 


2 a 7) = 
a= F : i cos() +c (n+ 5s) | are) 


where, A and nare the hardening parameters and @? is the equivalent plastic strain. However, this expression neglects the influence of 
the hydrostatic pressure and normalized Lode angle parameter on the hardening behaviour of materials. Therefore, Bai and co-workers 
proposed a generalized hardening rule with pressure and Lode angle dependence in the form: 


G=AZ [e + (l os) (set) 1) | (18) 


where cs constant defines the shape of the yield locus. For example, if the value of cz is equal to 1, the yield surface corresponds to von 
Mises yield criterion, while for cz = /3/2 it corresponds to Tresca yield criterion [40]. Thus, the modified Mohr-Coulomb model 


becomes: 
4/ : = cos() +¢ (1 + sin) | | (19) 


ae |: Coreen o)(seo 1) |. 


where c j, c 2 and cs are the three parameters of the modified Mohr-Coulomb model. This model was extended to enable an improved 
distinction between tension and compression stress states. In that case, the model considers four parameters: c}*, cj, c 1 and c 2, such 
that the equivalent strain to fracture is defined by: 


a F lore (es ~ 6) (seo) i) 


where the parameter c#* is defined as follows: 


i 1 for A>0 
Ge a 
A cy for 8 <0 


The parameter cj controls the asymmetry of the yield surface and consequently the fracture locus [71]; while c§ controls its 


1 


4/ : = cj cos) +c (1 + ssn) | | , (20) 


dependence on 6 [40]. To identify the parameters ci, and cj, the authors recommended the use of tests under the following conditions: 
axial-symmetric compression, plane strain and axial-symmetric tension. They mentioned that these parameters can also be determined 
from plasticity tests, since they are involved in the determination of the yield surface shape [40]. Nevertheless, the MMC model 
adopted in this work is the one given in Eq. (19), i-e., it has three parameters: c, c2, c3. Note that in this case, it is also necessary to 
identify the hardening parameters A, n, which can be determined using the stress-strain curve obtained from an uniaxial tensile test 


[2]. 
3.5. Lou and Huh, 2013 (L-H) model 


This fracture model describes the micro-mechanisms of nucleation of voids, growth and coalescence by the stress triaxiality, and the 
shear mode by the normalized maximum shear stress, Tmax, aS follows ([29,31)]): 


(=) (f + 3n) ya ae, (21) 
o 2 


where c 1, C 2 and c3 are model parameters and. 
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Tmax = (a = 63) (22) 
2 
Using the relationship between the principal stresses (01, 02,03), the Lode parameter, L (Eq. (8)) and stress triaxiality, 7 (Eq. (4)), the L- 
H, 2013 model can also be written as follows: 


() 


(tea) (Hest) #=> os 
since 

i _ —— (24) 
and 

Nl) =F = 14+ (25) 


where ( ) represents the Macaulay brackets. The authors initially suggested the introduction of a constant cut-off value for the stress 
triaxiality. However, they found that the cut-off value for the stress triaxiality can be affected by the temperature, the material 
microstructure ([72,31]) and the stress state, through the Lode parameter ([31,54]). Therefore, they proposed an extended version, 
defined as follows: 


Or) (esa) a. - 
with 
eer 
ALO =Z4C= nts aeatC (27) 


The two terms, 27,,,x/@ and f(y,L,C) are used to describe the shear voids coalescence and the voids growth, respectively. f(7,L,C) 
represents the dependency on L of the cut-off value for the 7. As previously mentioned, this cut-off value is also dependent on the 
material microstructure. Therefore, the C parameter was also introduced to define the height parameter of the cut-off value of 
triaxiality. This model encompasses a total of four parameters: c,, C2, cz and C. The fracture parameter C determines the ductile fracture 
occurrence for a certain loading condition. For instance, for C = 0, the cut-off value of the 7 is associated to the condition o, = 0 (see, 
Eq. (27)); while for C > 0, it corresponds to low stress triaxiality value. L-H model was successfully used in the numerical analysis for 
metal forming processes, under monotonic loading conditions, including shear dominated stress states, alike the MMC model [73]. 


3.6. Hosford-Coulomb, 2014 (H-C) model 


According to several experimental and theoretical studies, fracture model based on Mohr-Coulomb criterion combined with 
isotropic yield criterion enable accurate predictions of fracture under stress states ranging from shear to tension ([31,48,52]). 
Therefore, the M-C model was improved by substituting Tresca equivalent stress (see Eq. (13)) by Hosford’s one [47], to enhance the 
agreement between the experimental data and the model predictions [48]. The form of the M-C model expressed in terms of principal 
stresses is defined as follows: 


(0, — 03) +.¢1(0) + 03) = 2c2 (28) 


and the Hosford equivalent stress is given by 
ees é , ma 
Our = {5 (a — 62)° + (62 — 63)" + (0; — 03)" ) \ (29) 


where a presents the Hosford exponent. Note that the Hosford criterion becomes non-convex for a < 1, but there is no restriction with 
respect to convexity when it is used as localization criterion; thus, it is recommended to consider (0 < a<2) [48]. The H-C fracture 
model is defined in terms of the Hosford equivalent stress and normal stress as follows: 


Our + c(oy +03) = bo (30) 
= Cy _ __2¢9 
where c er and bo jut 
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As for the MMC model, this ductile fracture model can also be expressed as a function of the stress triaxiality, the Lode angle 
parameter and the mean stress, by using the geometrical construction of the principal stresses, deviatoric stresses and the Lode angle on 


the deviatoric plane, as shown in Fig. 1. 
The deviatoric principal stress can be expressed in terms of equivalent stress and Lode angle as follows: 


fs = [reos(o 

3 2. 2n 
on = i[22%0s (= _ 0) (31) 
Jo-vimels 

753 = jreo8(F ) 


Using the relationship between the principal stresses and the principal deviatoric stresses, oj = om + si, i = {1..3}, the principal 
stresses are defined as [40]: 


oO= (1 ee) on = (1 +A )on 
3n 1 


2. 
2e0s (= -0 f 
02 => 1+ On = (: + *) 
3n (32) 
2cos (— - 0) 
o3= 41+ dl Co= (1 +2) om 
3n n 


where either the Lode angle or the normalized Lode angle parameter-dependent functions [24] are used. The latter are defined as 
follows: 


2 1 
3° \6 
A) = 0s (= (3+ 8) ) where 9 = 1 —-— oy 


Then, by substituting the Eq. (32) into Eq. (30), the resulting Hosford-Coulomb fracture surface can be expressed in the (7, 0,6) space 
as: 


bo 
{1((i —f)" + hf)" + (i —f)") F + cn +f +f) 


a(n, 8) = (34) 


Fig. 1. Geometrical construction of the principal stresses, deviatoric stresses and the Lode angle on the deviatoric plane [26]. 
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The transformation of the Hosford-Coulomb criterion from the stress space (7,0,&) into the strain space (7,0, @;) requires the se- 
lection of an isotropic hardening law. If the simple isotropic power-law (see Eq. (16)) is considered, the H-C fracture model can be 
expressed as follows: 


nm 


4(1,8) = (%) (Ge A +A +B) ) +e +6) (35) 


Some authors proposed its generalization to allow the use of a more general isotropic hardening law, expressed as follows [47]: 


n 


(1,0) = (1 +0} ({u - fy b= Ay +h A) } HeOne +6) (36) 


In this form, the H-C ductile fracture model has three parameters: a, b and c, i.e. the hardening law parameters are identified in 
conjunction with the fracture ones, except for n which is assumed to be the equal to the one determined with the power isotropic law. 
This later form is the one adopted in this work. This model was calibrated using an extensive experimental database, including the 
analysis of localization and void fraction evolution. Besides, numerical simulations were performed to assess its predictive perfor- 
mances [74]. 


3.7. Quach et al., 2020 (Q-K) model 


Quach et al. claimed that most ductile fracture models cannot describe ductile fracture for various stress state, because they do not 
take into account all the macroscopic and microscopic mechanisms involved [9]. For instance, the empirical model proposed by Bao 
and Wierzbicki [53] considers the effect of the maximum principal stress, expressed in terms of the stress triaxiality, but neglects the 
effect of maximum shear stress, expressed by the Lode angle parameter. On the other hand, the MMC model does not include any 
micro-mechanisms. Therefore, Quach et al., proposed a new ductile fracture model (Q-K), taking into account the effect of both 
macroscopic and microscopic mechanisms on ductile fracture prediction, based on extensive studies and experimental observations. 
The Q-K model is expressed as follows ([9,58]): 


Cc 
@+9)°[G+ VIC) @) -G] 


o 


ef = 


(37) 


where Cj, C2 and Cs are the material parameters. This expression highlights that, from the macroscopic point of view, the proposed 
model considers the effect of both the maximum principal stress as well as the maximum shear stress. This fracture model is presented 
in function of the principal stresses, but using Eq. (32), it can be formulated as follows: 


C 


é:(n, L) = _———1—_ (38) 
aN) eat) eel) 
where 
3-L 1 
L)=n+ + 
ania 3V24+3 VP +3 
Rae Glee 


From the microscopic point of view, the influence of both maximum principal stress and normalized maximum shear stress on the 
void growth evolution are expressed in terms of both stress triaxiality and Lode parameter in the function g; |7,L]. The effect of sec- 
ondary voids and rotation of voids during void coalescence, which is governed by maximum shear stress, is considered as a function of 
Lode parameter, as established in the function g|L]. 


4. Experimental database and calibration procedure 


Two distinct materials were selected, namely 2024-T351 aluminium alloy and DP600 steel. For the 2024-T351, we revisit the 
experimental tests reported in [1,40,77,78] that were also broadly used to study the predictive ability of various ductile fracture 
models, in a wide range of engineering applications [35,8,9]. Therefore, this data could be well-established as a benchmark for 
calibration and verification of uncoupled ductile fracture models. For the DP600 dual phase steel, we revisit the experimental dataset 
obtained by Basaran [59], which was used to study the predictive ability of various ductile fracture models in the positive range of the 
stress triaxiality. In this study, the experimental data is applied also to compare the fracture predictive ability of seven uncoupled 
models for the same range of stress triaxiality. These materials were selected based on the diversity and abundancy of experimental 
results. 
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4.1. Experimental database: 2024-T351 aluminium alloy 


In order to cover a wide range of stress states, different samples were used for the 2024-T351 aluminium alloy, as shown in Fig. 2 
({2,35,79,80,54]). According to these latter references, the stress state changes during the test, although the test conditions attempt to 
impose monotonic loading paths. This presents a substantial issue for the calibration of the fracture models. For this reason, for each 
test, the stress state was analysed using a combined experimental and numerical approach [79]. Table 1 summarizes, for the 2024- 
T351, the average value of stress triaxiality and normalized Lode angle determined using this approach along with the equivalent 
strain to fracture. The set of 29 tests was divided into two subsets. The first subset includes tests performed under PSC (from #1 to #17) 
that was used for the identification of models’ ductile fracture parameters. The second subset comprises tests carried out in the positive 
stress triaxiality range (from #18 to 29), which was used to assess the fracture predictive ability of the calibrated models. 


4.2. Experimental database: DP600 steel 


The set of experimental tests carried out for the DP600 [59] used different types of samples, as shown in Fig. 3. Table 2 summarizes 
the average value of stress triaxiality and normalized Lode angle along with the equivalent strain to fracture. The set of 18 experi- 
mental tests was divided into two subsets: the first comprises tests performed under PSC (from #1 to #6) and was used for the cali- 
bration of the ductile fracture models; the second subset includes tests carried out in the positive stress triaxiality range (from #7 to 
#18) and was used to verify the accuracy of the calibrated models to predict failure. 


4.3. Calibration procedure of ductile fracture models 


The models were calibrated using a procedure to solve an optimization problem, in which the objective function corresponds to the 
gap between the experimental fracture strain and the value predicted by the ductile fracture model, expressed as follows: 


. 1 & (Zionainst — gait) 
ae eS 


(Wy Wa4--Wy) (Wy Was-Wn) i=I EFrxxpecimenial A 


(39) 


where (1,2, .--;Wn) denote the ductile fracture parameters, N is the total number of tests considered, and Efg,.cimenad ANA Efpredicread 
correspond to the experimental and predicted fracture strain values, for the same test i, respectively. 


Under plane stress condition 
(Tests #1 to #17) 
A 


Fig. 2. Geometries of the specimens used in the experimental tests for 2024-T351 aluminium alloy (adapted from [1,79]). 
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Table 1 

Experimental results of stress state and fracture strains for the 2024-T351 aluminium alloy ([1,40,77,78]). 
Test # Specimens n 7) L é ep 
1 Bi-axial compression —0.496 —0.398 0.367 —0.5852 0.349 
2 Cylinder compression (d/h = 0.5) —0.278 —0.8215 0.7946 —0.9609 0.4505 
3 Round notched, Compression —0.2476 —0.7141 0.6796 —0.9008 0.6217 
4 Cylinder compression (d/h = 0.8) —0.2339 —0.6809 0.6451 —0.8769 0.38 
5 Cylinder compression (d/h = 1.0) —0.2326 —0.6794 0.6435 —0.8758 0.3563 
6 Cylinder compression (d/h = 1.5) —0.2235 —0.6521 0.6155 —0.8543 0.341 
7 Pure Torsion 0 0 0 0 0.288 
8 Pure Shear 0.0124 0.0355 —0.0322 0.0557 0.2107 
9 Combined Shear and Tension 0.1173 0.3381 —0.3099 0.5064 0.2613 
10 Combined Tension and Torsion 0.1585 0.478 —0.4428 0.6822 0.2656 
11 Torsion with constant axial tension 0.191 0.556 —0.519 0.7664 0.253 
12 Plate with hole, Tension 0.3431 0.9661 —0.9594 0.9985 0.3099 
13 Pipe, Tension 0.3557 0.9286 —0.9155 0.9937 0.3255 
14 Dog-bone, Tension 0.357 0.9182 —0.9034 0.9917 0.4798 
15 Thin notched specimen (r = 4.763 mm), Tension 0.4319 0.696 —0.6607 0.8881 0.2441 
16 Thin notched specimen (r = 1.984 mm), Tension 0.4974 0.4486 —0.4145 0.6477 0.2204 
17 Thin notched specimen (r = 0.396 mm), Tension 0.5648 0.0942 —0.0855 0.1474 0.1951 
18 Solid square bar, Tension 0.3687 0.9992 —0.999 0.9999 0.3551 
19 Round bar, Tension 0.4014 0.9992 —0.999 0.9999 0.4687 
20 Notched bar(r = 14.503 mm), Tension 0.5243 1 —1 1 0.308 
21 Flat-Grooved, Tension 0.603 0.0754 —0.0684 0.1181 0.21 
22 Thick notched specimen(r = 13.970 mm), Tension 0.6052 0.1329 —0.1207 0.2072 0.1951 
23 Notched bar(r = 9.144 mm), Tension 0.6072 1 -1 1 0.2909 
24 Large notched bar (r = 12 mm), Tension 0.6264 0.9992 —0.999 0.9999 0.283 
25 Notched bar(r = 6.096 mm), Tension 0.6765 al -1 1 0.2727 
26 Thick notched specimen(r = 4.763 mm), Tension 0.7304 0.0692 —0.0062 0.1084 0.2017 
27 Notched bar(r = 4.470 mm), Tension 0.7631 1 —1 1 0.2447 
28 Notched bar(r = 3.048 mm), Tension 0.8556 1 —1 1 0.2022 
29 Small notched bar (r = 4 mm), Tension 0.9274 0.9984 —0.998 0.9999 0.1665 


The objective function is minimized using the generalized reduced gradient (GRG) method for solving nonlinear problems, 
implemented in Excel. This method uses the gradient of the objective function as input to change the models’ parameters until it attains 
a minimal value for the partial derivatives, i.e., an optimal solution is attained. The advantage of GRG nonlinear solver is its fastness; 
however, the main drawback is the dependence on the initial solution and the menace of falling into a local minimum, close to the 
initial input, instead of converging into the global optimal solution. To overcome this issue, the GRG multistart option was activated, in 
order to let the algorithm create a randomly distributed population of initial input values. By starting several times from different 
initial inputs, there is a much higher probability that the obtained solution falls into the global optimum. 

The choice of the this calibration procedure is based on the fact that it enables achieving good results, when it was firstly used by 
Bao et al. [1] to calibrate fracture models for the 2024-T351 aluminium alloy and, subsequently, by Bai and Wierzbicki [35] as well as 
other researchers (see [8,81,82,23,9,45,83,46]). Thus, the performance of this calibration approach was satisfactorily demonstrated 
for distinct fracture models. 

As previously stated in sections 4.1 and 4.2, the set of the experimental tests listed in Table 1 and 2, was split into two subsets. For 
2024-T351 aluminium alloy, the first subset corresponds to 17 experimental tests (from #1 to #17), allowing to cover a range for 7 
from —0.496 to 0.5648, and for € between —0.9609 and 0.9985. For DP600 steel, the first subset contains 6 tests (from #1 to #6) 
corresponding to the range from —0.06 to 0.66 for 7 and to the range from —1 to 0.96 for &. According to the PSC, the relationship 
between y and @ is established by Eq. (9) or between 77 and é by Eq. (10) [2]. The use of his relationship simplifies the analysis of the 
fracture models. Thus, the fracture locus can be represented in a 2D space:(n, é;). 

Depending on the model considered, the number of damage parameters varies from three to six, as listed in Table 3. However, some 
models also require the identification of the hardening law parameters (e.g., X-W, MMC and H-C models). In this work, a power 
hardening law was considered as expressed by Eq. (16)). The hardening parameters were identified using the stress-strain curve of the 
tensile test performed with the round bar and with the un-notched flat specimen, for 2024-T351 and DP600 steel, respectively ([2,59]). 
The obtained hardening parameters for 2024-T351 aluminium alloy and DP600 steel are A = 740 MPa, n = 0.153, and A = 918.87 
MPa, n = 0.16, respectively. For the X-W model, the even integer closest to 1/n, corresponding to m, is taken equal to 6, for both the 
2024-T351 and the DP600. 

It is noteworthy that the L-H model is 4 parameter-depend. The fourth parameter is denoted by C and it can also be identified using 
the optimization procedure [54]; however, it is recommended to choose the parameter C such that the 77 value for all fracture tests is 
higher than the cut-off value for 7 predicted by Eq. (27) [54]. In order to discuss the selection of this parameter, the influence of C on 
the cut-off value is analysed taking into account Eq. (27). Fig. 4 depicts the evolution of the cut-off value of 7 as a function of the Lode 
parameter considering different values of C, determined when Eq. (27) is equal to zero. Fig. 4 shows that the cut-off value of 7 decreases 
with the increase of C, such that when the value of Cis higher than 1/3 (e.g., C = 0.4 or C =0.5), the cut-off value of 7 corresponding to 
Lode parameter L = 1 will be lower than —2/3. However, all tests considered are under PSC, which implies that 7 is only varying 
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Fig. 3. Geometries of the specimens used in the experimental tests for the DP600 steel ([59]). 


Table 2 
Experimental results of stress state and fracture strains for the DP600 steel ([59]). 


Test # Specimens UT 0 L é e 


&f 

1 Butterfly —10°, combined compression-shear —0.06 —0.1608 0.146 —0.25 1.02 
2 Butterfly 0°, shear 0.07 0.1740 —0.158 0.27 1.08 
3 Butterfly 10°, combined shear-tension 0.18 0.4936 —0.458 0.7 0.7 

4 Butterfly 20°, combined shear-tesnion 0.34 0.8193 —0.792 0.96 0.95 
5 Nakazima 1 0.64 —0.5303 0.494 —0.74 0.55 
6 Nakazima 2 0.66 -1 1 -1 0.84 
7 Smooth flat, tension 0.61 0.5209 —0.485 0.73 1.41 
8 Notched flat 1 R = 4 mm, tension 0.63 0.4017 —0.37 0.59 0.81 
9 Notched flat 2 R = 2 mm, tension 0.67 0.2344 —0.214 0.36 0.78 
10 Notched round R = 4 mm, tenion 0.7 1 -1 1 1.19 
11 Notched round R = 2 mm, tension 0.74 1 -1 1 1.08 
12 Grooved flat R = 4 mm, tension 0.76 0.0063 —0.006 0.01 0.98 
13 Grooved flat R = 2 mm, tension 0.8 0.0063 —0.006 0.01 0.84 
14 Butterfly 60°, combined shear-tension 0.82 0.1806 —0.164 0.28 0.44 
15 Notched round R = 1 mm, tension 0.82 1 -1 iT 0.88 
16 Grooved flat R = 1 mm, tension 0.88 0.0063 —0.006 0.01 0.7 

17 Notched round R = 0.5 mm, tension 0.95 1 -1 1 0.67 
18 Grooved flat R = 0.5 mm, tension 1.03 0.0063 —0.006 0.01 0.59 


between —2/3 and 2/3 and the Lode parameter L is in the interval [-1, 1]. 

For the 2024-T351 aluminium alloy, according to test #1 (see, Table 1), the cut-off value of 7 should be lower than —0.496, forL = 
0.367. This only justifies the adoption of the C parameter slightly higher than zero, to be certain that the fracture in compression 
conditions will occur for this value. Due to the lack of data, it was decided to consider C equal to 1/3 to assure that the cut-off value for 
the 7 will be lower or equal to —2/3. Unfortunately, for the DP600 steel, there is no experimental data for compression conditions. 
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Table 3 
Damage parameters considered in the calibration procedure for each model. 
Models Number of damage parameters Damage parameters 
B-W 6 D,,D2,D3,D4,D5,D6 
M-W 4 D,,D2,D3,D4 
X-W 4 C1,C2,C3,C4 
L-H 3 €1,€2,C3 
MMC 3 C1,C2,C3 
H-C 3 a,b,c 
QK 3 Ci,C2,C3 


1.0 -0.8333 -0.733 -0.666 -0.533 -0.433  -0.333 


o,= 0 for C=0 


0.0 


Lode parameter 


-0.5 


Cut-off plane of stress triaxiality 


i 3-L 
S(n.b,C)=2+C =n +a +t C=0 
Cc 3V ; 


°1.2 -1.0 -0.8 -0.6 -0.4 ~0.2 0.0 


Stress Triaxiality 77 


Fig. 4. Cut-off plane of 1 for different values of C for the Lode parameter L between —1 and 1. 


Thus, the same value was assumed for C parameter to ensure the targeted consistency for the calibration strategy adopted. 


4.4. Choice of stress states for models calibration and verification 


The main purpose of this study is to assess the effect of using different type and number of tests on the identified models’ parameters 
and then to evaluate the predictive ability of these models. In this sense, the available experimental data, for the 2024-T351 aluminium 
alloy, was split in different groups. Note that, the minimum number of experimental tests to be used for the calibration of the more 
flexible fracture model is 6 (see Table 3). Thus, groups of 6 tests were selected for which the stress states covered were carefully chosen 
in each group. Then, for comparison purpose with the calibration results obtained using the different groups of tests, the whole 
experimental dataset (Group 4) was used in the calibration of the models as routinely performed in previous works. 

The different sets of 6 tests were carefully elicited from the 17 experimental tests available. The first calibration strategy uses the 6 
following tests: #7, #8, #11, #14, #16, and #17, for which 7 ranges from 0 to 0.5648 and é is between 0 and 1. This is referred as 
Group 1. The second is referred as Group 2 and it is composed of the following 6 tests: #2, #6, #8, #10, #13 and #16, as listed in 
Table 1, for which 7 ranges from —0.278 to 0.4974, and é is between —0.9609 and 0.9937. These 6 tests were selected to analyse the 
influence of a wider amplitude foré. The third is referred as Group 3 and includes the following 6 experimental tests: #3, #4, #7, #11, 
#14 and #17, as shown in Table 1. In this group, the 7 ranges from —0.2476 to 0.5648 and é ranges from —0.908 to 0.9917. Therefore, 
Group 2 covers a slightly lower range of 7 when compared to Group 3 and a similar range for €. More importantly, they present 
dissimilar equivalent strains to fracture, particularly between tests (#2, #3) and (#13, #14). Therefore, it allows to study the influence 
of considering different types of tests, for the same range of stress states. The last group, referred as Group 4, uses the full set of 17 
experimental tests. This means that the whole range of available 7 and é values is covered. Thus, this group shows a slight increase of 
the 7 range, when compared with the ones used in the previous groups and a similar range for é. 

For the DP600 steel, due to the lack of experimental data, only one group of 6 tests is considered in the calibration procedure. The 
group of 6 tests are labelled from #1 to #6 (see, Table 2). This means that they cover a range for 7 from —0.06 to 0.666 and for € from 
—1 to 0.96. The comparison with the ranges considered in the groups used for 2024-T351 show that the more similar is Group 2, 
although it covers a slightly higher range of 7. Group 1 presents a more similar range for 7, but a lower range for &. Therefore, both 
groups are analysed in the comparative study concerning the predictive ability of the fracture models, for both materials. 

The predictive ability of the calibrated fracture models was assessed based on the comparison between predicted fracture strain and 
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the observed one. For the 2024-T351 aluminium alloy, the subset of tests from #18 to #29 is considered. For the DP600 steel, the 
subset of tests from #7 to #18 is applied. For both materials, these subsets correspond to tests performed in the high 7 range, i.e., 
n21/3. According to many authors (see e.g. [40,31]), the fracture occurrence in smooth and notched round bars or flat specimens, 
subjected to tension, is governed by the void growth and linkage of micro-cracks mechanisms. However, in compression, the shear 
mode dominates. Unfortunately, such fracture mechanisms are not properly captured by most of the ductile fracture models. In this 
context, it is important to assess the accuracy of the fracture prediction under the upper bound assumption, i.e., = +1 and in the 
triaxial loading state, as occurs for the experimental set of tests from #18 to #29, for 2024-T351, and from #7 to #18 for DP600. Note 
that the test subjected to bi-axial compression loading (test #1) was considered neither in the calibration nor in the verification of the 
prediction ability of M-W model, because of the existence of a cut-off value of 7 equal to—1/3. 


5. Results and discussion 


This section starts with the analysis of the influence of the set of data selected on the calibration and predictive ability of the seven 
uncoupled models. This study is performed for the 2024-T35 aluminium alloy. Afterwards, the results obtained for both materials are 
analysed in order to assess to what extent the conclusions drawn for 2024-T351 can be extended to other materials. The average error 
values presented, either for the calibration or for the verification of the fracture models prediction ability is based on Eq. (39). 


5.1. Analysis of the effect of the experimental dataset on the fracture models for 2024-T351 


5.1.1. Analysis of the calibration error 

As previously mentioned, four groups of tests were used to calibrate the seven fracture models, Group I, where I = {1, 2, 3, 4}. 
The identified damage parameters are listed in Table 4. Groups 1, 2 and 3 consist of only 6 tests, while Group 4 considers 17 tests in the 
calibration procedure. Therefore, the average calibration error is calculated considering different types and number of tests. For all 
Groups. the calibration error is obtained at end of the optimization procedure, i.e. corresponds to the minimum of the objective 
function (Eq. (39)) for the specific set of tests considered. Nevertheless, for Groups 1, 2 and 3, the whole set of 17 tests is also used in the 
calculation of the average error. In other words, 11 tests among the 17 tests, which are not used in the calibration procedure, are 
considered. This allows to evaluate the impact of neglecting specific types of tests in the calibration procedure. 

In Fig. 5, it is observed that the average calibration errors obtained for Group 1 and 2 are relatively lower than the ones obtained for 
Group 3. In addition, for all models, the average calibration errors for Groups 1, 2 and 4 are lower than 9%, except the M-W. Besides, 
for all models, the Group 2 has lower average calibration errors than Group 1, except for the M-W and the Q-K models. 

It is somehow expected that the average error calculated considering tests not used in the calibration procedure is higher than the 
objective function value at the end of the optimization procedure. The magnitude of this error depends on the Group used for the 
models’ calibration, as depicted in Fig. 5. It is interesting to note that, from the 6 tests sets, Group 2 leads to the lowest values for the 


Table 4 
Damage models parameters identified for 2024-T351 using Group 1, Group 2, Group 3 and Group 4 of experimental tests. 

Model Damage Parameters 
Parameters Group 1 Group 2 Group 3 Group 4 

B-W Dy 0.5451 0.5180 0.5242 0.4580 
Dz 0.001 1.1550 0.0001 1.0857 
D3 0.21 0.2091 0.2875 0.2196 
D4 0.0707 0.1342 0.7464 0.2984 
Ds 1.1476 0.4264 0.4426 0.342 
De 0.6333 1.0924 0.0001 1.7676 

M-W Dy 0.2069 0.1123 0.1597 0.1135 
Dz 0.6070 1.2256 0.3224 1.2259 
D3 —2.937 —4.99 —2.171 —5.0 
D4 0.079 0.1185 0.0 0.1185 

X-W Cy 2.8961 0.7227 0.9151 0.7062 
C2 3.1124 1.2151 0.0001 1.1565 
C3 0.2110 0.2677 0.2817 0.2732 
C4 0.1392 0.4101 0.6505 0.4599 

L-H C=1/3 Cy 7.7668 3.947 7.3072 3.8076 
C2 0.2984 0.3127 0.1810 0.3309 
c3 0.5773 0.3591 0.5711 0.3540 

MMC cy 0.0057 0.0288 0.0174 0.0377 
c2(MPa) 339.849 324.95 345.84 326.9 
c3 1.0168 0.9477 1.0227 0.9467 

H-C a 0.9285 1.1621 1.4942 1.2351 
b 0.5897 0.3492 0.2596 0.334 
c 0.0137 0.0303 0.0461 0.0352 

Q-K Cy 0.3770 0.4711 0.401 0.4406 
C2 0.2107 0.5836 0.3318 0.4348 
C3 7.2130 4.1789 7.1962 3.5115 
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average error calculated with the 17 tests, which is also quite similar to the one calculated for Group 4. One could relate this to the wide 
range for the variables that define the stress states of the tests. However, group 2 and 3 presents a quite similar range for the stress state 
variables. 

In order to better understand these results, Fig. 6 displays the relative error associated to each experimental test used in the 
calibration. Fig. 6 (a) shows for Group 1, that all models fail to describe the fracture strain for the pure torsion test (test #7), for which 
the relative errors are higher than 25%. Nevertheless, all models predict the fracture strains in the pure shear test (test #8) with an 
error lower than 10%. In this context, it is interesting to mention that these two experimental data points, corresponding to pure 
torsion (test #7) and pure shear (test #8) tests, present similar values for the stress state variables; however, the equivalent plastic 
strain values are slightly different (7.7% difference between the highest and the lowest). 

In fact, the comparison with the relative errors obtained with Group 2, shown in Fig. 6 (b), indicates that it is the inclusion of test #7 
that is contributing to the higher average calibration errors obtained with Group 1. These results seem to highlight the negative impact 
of the inclusion of tests with similar stress triaxiality and Lode parameter, but different fracture strains, which induce doubts about the 
fracture strain. However, it is noteworthy to remind that the measure of fracture strains is marred by uncertainty related to numerous 
sources. Fig. 6 (b) also shows that when using Group 2, there are no particular difficulties in capturing any particular test. 

Fig. 6 (c) shows that, when using Group 3 for the models’ calibration, all models, except M-W, present difficulties in capturing the 
compression of a round notched specimen (test #3), for which the relative errors are close to 40%. This indicates that it is the inclusion 
of this test that is contributing to the higher average calibration errors obtained with Group 3 (see Fig. 5). In this context, it is 
interesting to note that Group 4 also presents relative errors close to 40%, for all models. The fact that the average calibration error is 
lower, when compared with the one obtained with Group 3, is a consequence of the total number of tests considered (see Eq. (39)). 

In order to further elucidate on the contribution of each test to the average error calculated with 17 tests, Fig. 7 presents the relative 
error values obtained for each test, performed under PSC, which has not been considered in the models’ calibration for Group 1, 2 and 
3. It is noteworthy that Group 1 globally presents high relative errors for tests with € < 0, which is probably related to the range 
considered for this parameter in the calibration set. Nevertheless, it also leads to worst predictions when comparing with Group 2, for 
tests with € = 1. In fact, Groups 1 and 3 lead to similar predictions for tests #12 and #13 (Plate with hole and pipe, under tension). It is 
also interesting to note that Groups 1 and 3 lead to similar average error values (calculated with all the tests in Fig. 5), for L-H, MMC, H- 
C, and Q-K models. However, the same average error is higher for models that are slightly more flexible, i.e., B-W, M-W and X-W 
models. For these models, the average errors obtained with the Group 1, when considering all tests under PSC are higher than 35%, as 
illustrated in Fig. 5. As previously mentioned, this high average error can be explained by the absence of experimental tests with € < 0 
in the calibration procedure (tests from #2 to #6, in Fig. 7 (a)). Therefore, B-W and X-W models seem to present a higher sensitivity to 
the stress state range considered in the calibration than the other models. 

The sensitivity of the more flexible models to the dataset selected for calibration is particularly evident for the B-W, which attains 
an average calibration error close to zero only when Group 2 is selected. Nevertheless, as for Group 4, all models are not able to predict 
failure in the round notched compression test (test #3) and the tension test carried out on the dog-bone specimen (test #14). Fig. 5 
highlights that the M-W model always leads to the highest values of the average calibration errors. 

The analysis of Fig. 6 indicates that this model always presents more difficulties in capturing some of the calibration test values. For 
instance, when using Group 2, the relative errors values predicted by M-W model for the tests undergoing compression loading (test 
#2) and pure shear (test #8) are higher than 40% (Fig. 6 (b)), which obviously contribute to the increase of the average error value, 
when compared with other models. Note that for M-W model, the existence of cut-off value at the negative range of stress triaxiality 
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Fig. 5. Average errors obtained in the calibration procedure for 2024-T351: considering 6 tests in case of Group I, I= {1, 2, 3} and 17 tests for 
all groups. 
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Fig. 6. Relative errors obtained in the calibration procedure performed for 2024-T351using 6 tests: (a) Group 1; (b) Group 2; and (c) Group 3; and 
17 tests (d) Group 4. 


and its poor flexibility (only one damage parameter), limit its predictive capacity in this stress state range. Thus, when using Group 3, 
the M-W model fails to predict the fracture in the compression test on short cylinder specimen (test #4), pure torsion test (test #7), but 
also the tensile test carried out on the dog-bone specimen (test #14). 

Fig. 8 shows the 2D fracture loci associated to the seven models that were calibrated using Group I, withI = {1, 2, 3, 4}. In this 
figure, the tests used in the calibration procedure, for each Group, are marked with white-coloured diamonds. Globally, for the same 
calibration group, the trend of the fracture strains obtained with the seven selected models are very similar. These is particularly 
evident for L-H, MMC and Q-K models. In fact, the group selected seems to have a bigger impact in the local peaks for 7 = —1/3 andy = 
+1/3 (see zone A and zone B in Fig. 8, respectively). Fig. 8 (a) and (b) show a significant difference in the fracture locus obtained using 
Group 1 and 2, for all models, especially in zone A and zone B. This is certainly related with the fact that Group 1 considers only the 
positive range of 7 and Lode angle parameter. Significant differences are also observed in the 2D fracture locus obtained by the selected 
models when considering the two groups with very similar ranges for the stress states but slightly different values of fracture strain, i.e. 
Group 2 and 3. The three tests (tests #3, #7, and #14 used by Group 3) present clearly higher values of equivalent strain to fracture, 
when compared with the neighbouring stress states (tests #2, #8, and #13 used by Group 2). Fig. 8 (c) shows that Group 3 leads to a 
shift of the 2D fracture locus upward, with a strong change in its shape, for all models, except B-W and H-C. The B-W, X-W and H-C 
models present a quite different trend for 7 =—1/3 and y = +1/3 (see zone A and B in Fig. 8 (c)). For the B-W model, the calibration 
using Group 3 shifts the fracture locus downward in the negative range of 7. However, it presents a similar evolution of the equivalent 
strain to fracture as the other models, in the high range of 7. Furthermore, the X-W model predicts higher values of the fracture strains 
for 7 =—1/3 and 7 = + 1/3, whatever the calibration group adopted. Moreover, it attains the highest magnitude for the fracture 
strain for 7 =+1/3 when using Group 1 (zone B in Fig. 8 (a)). Thus, this model shows a higher sensitivity to the set of tests selected for 
the calibration than the other models. The fracture loci predicted by the H-C model using Group 3 is quite different of the ones obtained 
with the other models, since there is no strong increase of the fracture strain in the low range of 7 (0<7<1/3), as shown in Fig. 8 (c). It is 
clear from Fig. 8 (b) and (d) that Group 2 and Group 4, lead to a similar trend of 2D fracture locus, except for B-W and Q-K models. The 
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Fig. 7. Relative errors obtained in the calibration procedure performed for 2024-T351using 6 tests, for the remaining 11 tests performed under PSC: 
(a) Group 1; (b) Group 2; and (c) Group 3. 


fracture locus predicted by the B-W model with Group 2 is slightly concave for high 7 values while it becomes a little convex when 
using Group 4, particularly for 7>1/3. The Q-K model presents a quite different trend of the fracture locus, for the negative range of 7 
(see zone A in Fig. 8), when comparing Group 2 with Group 4. The calibration using the 6 tests (Group 2 in Fig. 8 (b)) shifts the fracture 
locus upwards, for 7 = —1/3. The shift downwards of the 2D fracture locus observed for Group 4 seems to be related with the larger 
number of tests close to negative stress state values (tests #1 and #3-5). 

In summary, the same calibration procedure performed with either different types and/or numbers of tests can lead to different 
average calibration errors, which tends to be lowest for the models that are more flexible. These models are also more sensitive to the 
range considered for &. Significant differences are also observed in the fracture locus obtained using different sets in the calibration. In 
this context, the fracture locus shape is also affected by the type of tests selected, within a similar range for 7 and é. 


5.1.2. Analysis of the prediction ability 

In this section, models’ prediction ability is assessed based on the comparison of the relative errors between the predicted fracture 
strain and the one experimentally measured, for the 12 tests from #18 to 29#, carried out in the high 7 range (not considered for the 
models’ calibration). Fig. 9 shows the relative error calculated for each test, considering the seven ductile fracture models calibrated 
using the different calibration sets. 

Globally, Group 2 and 4 lead to similar distributions of the relative prediction error, which could be expected based on the analysis 
performed in the previous section. Groups 1 and 3 lead to similar levels of prediction accuracy, except for the X-W and H-C models, as 
shown in Fig. 9 (a) and (c), respectively. In fact, it is interesting to note that the M-W model is the only one presenting relative 
prediction errors lower than 50%, whatever the set adopted in the calibration procedure. Nevertheless, this is certainly related with its 
lower flexibility. On the other hand, the H-C model is the one presenting the higher sensitivity in the relative prediction error to the set 
used in the calibration procedure. In fact, when using Group 3, the H-C model is the one presenting the best prediction, followed by M- 
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Fig. 8. Comparison of 2D fracture locus obtained with the seven uncoupled models for the 2024-T351: (a) Group 1; (b) Group 2; (c) Group 3; and 
(d) Group 4. 


W model. All the other models fail to predict the fracture for the upper bound assumptions, i.e., for € = + 1, when using Group 3. In 
this context, it is interesting to note that the type of tests seems to have a greater influence on the predictive ability than the range of 
stress states considered. In order to better understand the differences in the predictive ability of each model, when considering 
different sets in the calibration, the 3D fracture loci (plotted in the 3D space (1, £,@s)) are also presented. Nevertheless, for each model, 
only the 3D fracture loci that lead to clearly distinct relative error predictions are shown, to simplify the discussion. In these figures, the 
2D fracture locus is presented by the thick dash line, while the normalized uniaxial compression, plane strain, and uniaxial tension 
conditions are represented by the red, blue and black-coloured solid lines, respectively. The experimental tests used in the calibration 
procedure are marked with white-coloured diamonds. The tests performed under PSC, that were not included in the calibration 
procedure are marked by black-coloured circles. The tests used for the prediction accuracy, from #18 to #29, are displayed with the 
red-coloured circle marks. The fracture strain predicted for each test condition (7, €) defined in Table 1 is marked with blue-coloured 
circles. The relative errors for each test (#1 to #29) are denoted by black thin lines. 

Regarding the B-W model, as previously mentioned, Groups 2 and 4 lead to quite distinct 2D fracture loci, particularly for the high 
range ofy. Therefore, the shape of the 3D fracture loci is quite similar (concave shape), for low 7 values; but it becomes progressively 
flattened, for 2024-T351 Group 2, and a little convex, for 2024-T351 Group 4, at high 7 values (see zone C in Fig. 10 (b) and (c), 
respectively). The 3D fracture locus obtained with Group 3 has a concave shape for a wide range of 7 (Fig. 10 (c)). However, it is 
flattening out, for negative range of 7. The 3D fracture locus is always concave when considering Group 1 (see Fig. 10 (a), respectively). 
It should be also mentioned that the calibration set also has a strong influence on the magnitude of the equivalent fracture strain, 
particularity in the lower and the upper bounds, i.e. for € =—1 and é = +1. 

Although the four calibration sets are based on distinct conditions, the shape of the fracture surface predicted by M-W model is 
quite similar. Therefore, Fig. 11 presents only the one obtained with Group 4. The fact that M-W model is independent and the 
existence of the cut-off value at the negative 7 limits its prediction capability (see also Fig. 9 (b)). Nevertheless, the magnitude of the 
fracture strains is slightly affected by the calibration set, as shown in Fig. 8. 

Fig. 12 reveals some differences in the 3D locus shapes obtained with different calibration sets, for the X-W model. Note that groups 
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Fig. 9. Relative errors of fracture models prediction for the 12 tests #17- #29 using (a) 2024-T351 Group 1, (b) 2024-T351 Group 2, (c) 2024-T351 
Group 3, and (d) 2024-T351 Group 4. 


2 and 4 present similar trends of the 2D fracture loci and, consequently, also identical levels of prediction accuracy. Whatever the 
calibration set adopted, the 3D fracture loci of the X-W model has a concave shape for low 7 values. For high values, the shape changes; 
when Group 3 is adopted it remains concave, while it flattens out when using Group 4, and becomes convex when considering Group 1, 
(see zone C in Fig. 12). It is also interesting to note that, as for the B-W model, the calibration set also has a strong influence on the 
magnitude of the equivalent fracture strain, particularity in the lower and the upper bounds, i.e. foré =—land&= +1. 

Fig. 13 presents the 3D locus obtained for the H-C model, when using Group 1, 3 and 4. In this case, the bigger differences are 
observed when Group 1 is considered; i.e. the calibration procedure considers only tests in the positive range of 7. As for B-W and X-W 
models, Group 1 leads to a concave fracture locus. The change of the calibration set results in the flattening out of the fracture locus, 
including for € = 0, and to a decrease of the fracture strain, particularly for € = +1 (see Fig. 13). Thus, the 3D fracture locus of the X-W, 
B-W and H-C models show high sensitivity to the type and number of tests used in their calibration, which leads to different prediction 
errors, as shown in Fig. 9. The comparative study of the MMC, L-H, H-C and Q-K models shows a similar level of prediction accuracy 
under different loading conditions, when considering different number of tests in their calibration procedure, i.e., Group 2 and 4, as 
previously shown in Fig. 5 and Fig. 9. These models, except H-C, also show quite similar results when adopting Group 1 and 3, as 
depicted by the average errors obtained with the 17 experimental tests in Fig. 5. Therefore, for the MMC, L-H and Q-K models the 3D 
fracture loci are presented only for Group 3 and 4, in Figs. 14, 15 and 16, respectively. The analysis of these figures shows that the 3D 
surface of these models is quite similar when the calibration is performed with the 6 experimental tests, considered in Group 3; it 
presents always a concave shape. For MMC, L-H and Q-K models, the usage of a large number of tests (17 tests) leads to a shift of the 
boundary curves corresponding to € = —1 and € = +1 downwards, as shown in Figs. 14, 15 and 16, respectively. Moreover, their global 
shape is flatter. 

The analysis of region D (see Fig. 14 (b)), confirms that the B-W, MMC, and H-C models show the more accurate predictions at € = 

+1, when using Group 4. In contrast, L-H and Q-K models have higher average prediction errors in the high 7 range. In fact, in this 
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Fig. 10. 3D fracture loci obtained for the B-W model considering different sets of tests in the calibration procedure performed for the 2024-T351 
aluminium alloy: (a) Group 1; (b) Group 3; (c) Group 3 and (d) Group 4. 


range, it is found that the effect of £ on the 3D fracture loci predicted by these models becomes very weak, especially for 720.5 (see 
Zone C, as depicted in Fig. 10 (c)). 

In summary, the comparative study on ductile fracture prediction for the 2024-T351 aluminium alloy shows that the B-W model 
provides the best prediction accuracy, followed by X-W, MMC and H-C models. These models show a very similar level of prediction 
accuracy when using sets with different number of tests, i.e., Group 2 and 4. Moreover, L-H and Q-K models show satisfactory results in 
the negative and low-v ranges, but they fail to predict the fracture in normalized uniaxial tension. Other studies have shown that Q-K 
model is unable to predict failure under tension, although the authors state that its development is based on the mechanisms of 
nucleation, void growth and coalescence ([31,72]). 

The results obtained with the calibration performed using Group 1, which includes only tests in the positive 7 range, exhibits a 
similar level of prediction for X-W and M-W models, followed by B-W, MMC, L-H and H-C models. When using Group 3 in the cali- 
bration the H-C model shows the best prediction, followed by M-W model. It can be concluded that the different types and number of 
tests used in the models’ calibration have an influence on the shape of the 3D fracture locus, particularity its boundary curves, cor- 
responding to € 1,é=Oandé = +1. Consequently, they have a strong impact on the prediction accuracy of models, except for 
the M-W model, due to its poor flexibility. 
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Fig. 11. 3D fracture loci obtained for the M-W model for the 2024-T351 aluminium alloy when using Group 4. 
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Fig. 12. 3D fracture loci obtained for the X-W model considering different sets of tests in the calibration procedure performed for the 2024-T351 
aluminium alloy: (a) Group 1; (b) Group 3and (c) Group 4. 


5.1.3. Comparative study with previous works 
Some of the fracture models selected have been calibrated and evaluated in previous studies, also using the experimental data 
available for the 2024-T351. In this context, it is interesting to shed some lights on the similarity and the fundamental disagreement 
between the current work and the previous works, where the calibration procedures used the, as will be shown in the following. 
The parameters of the MMC, H-C and L-H models were calibrated, using the least square method, applied to the same 17 exper- 
imental tests used in Group 4. The trends of the fracture locus obtained as well as the average calibration errors are very similar to the 
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Fig. 13. 3D fracture loci obtained for the H-C model considering different sets of tests in the calibration procedure performed for the 2024-T351 
aluminium alloy: (a) Group 1; (b) Group 3and (c) Group 4. 


ones reported in this study. Nevertheless, the predictive ability is quite different for tests #19, #20, #22, #24 and #28; it is better in 
the current study for tests #20, #22 and #24 and worst for tests #19 and #28, which can be related with the different methodology 
adopted to solve the non-linear problem [8]. The Q-K and MMC models were also calibrated and compared with the L-H one (see 
details in [54]) by Quach et al. [9]. In this case, a set of 18 tests was used in the calibration procedure, by adding to the 17 tests, 
considered in Group 4, the tensile test on a flat grooved specimen (test #21) (see Table 1). Note that Group 4 leads to a good prediction 
of this test, which is not considered in the models’ calibration (see Fig. 9 (d)). The results obtained with 18 tests present a shift of the 2D 
fracture locus upwards, particularity for the local peaks corresponding to ny = —1/3 and 7 = + 1/3, when compared with the results 
obtained for Group 4; in both studies the results for the three uncoupled models are quite similar. This shift corresponds to a lower 
relative error for tests #3 and #14 than the one obtained in this work with Group 4. By contrast, the relative errors of the tests #2, #4, 
#5, #6, #12 and #13 predicted with Group 4 are lower, for L-H and MMC models. For Q-K model, only the relative errors corre- 
sponding to tests #12 and #13 obtained with Group 4 are lower than the ones reported in [9]. The shape of 3D fracture loci is similar 
for normalized uniaxial tension, but it becomes different for normalized uniaxial compression, which corroborates the impact of the 
method selected to solve the non-linear problem. In another study, the X-W and MMC models were calibrated using 12 tests, corre- 
sponding to the range from —0.278 to 0.603 for 7 and to the range from —0.969 to 0.99 for € [26]. Note that this set includes two tests 
(#18 and #21) that were not included in the current study, while it discards tests #1, #7, #10, #11, #15, #16 and #17 (see Table 1). 
Compared with Group 4, the set of 12 tests has a similar range for € and a slightly lower range for 7. The damage parameters and the 
shape of 3D fracture loci obtained with this set for MMC model are similar to the ones obtained with Group 4. For the X-W model, 
although the shape of the 3D fracture locus obtained with 12 tests and Group 4 is quite similar for the negative range of 7, it becomes 
different (convex shape for the calibration with 12 tests, as occurs in the current study for Group 1 (see Fig. 10 (a)), for positive high 7 
values ([2,26]). These results confirm that more flexible models are more sensitive to the type and number of tests considered in 
calibration. Finally, the B-W model was also calibrated using the same type of tests using a methodology based on the separate 
identification of the damage parameters for each of the three boundary curves, i.e., € = —1, =O and é = +1 ([32,35]). The set used 
consists of 15 tests, ten coincident with the ones used in Group 4 (tests #2 - #6, #8, #9 and #12 - #14) and other 4 that were not (tests 
#18, #19, #21, #24 and #29). This corresponds to a range for 7 from —0.278 to 0.93 and for € between —0.91 and 1. As results, the 
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Fig. 14. 3D fracture loci obtained for the L-H model considering different sets of tests in the calibration procedure performed for the 2024-T351 
aluminium alloy: (a) Group 3 and (b) Group 4. 
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Fig. 15. 3D fracture loci obtained for the MMC model considering different sets of tests in the calibration procedure performed for the 2024-T351 
aluminium alloy: (a) Group 3 and (b) Group 4. 


shape of 3D fracture surface obtained with 15 tests in the calibration seems very similar to the 3D fracture locus obtained with Group 2 
and 4, for =Oand é = +1, butit becomes a slightly different at € = —1, for the negative values of stress triaxiality. This shows that 
not only the number and type of tests has an impact on the calibration of B-W model but also the methodology considered in the 
optimization process. 
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2024-T351 Group 4 Q-K 


Fig. 16. 3D fracture loci obtained for the Q-K model considering different sets of tests in the calibration procedure performed for the 2024-T351 


aluminium alloy: (a) Group 3 and (b) Group 4. 


5.2. Analysis of the fracture models for DP600 steel 


5.2.1. Analysis of the calibration error 


Table 5 lists the ductile fracture parameters obtained for the seven models selected, for DP600. Fig. 17 shows the percentage 
average calibration error calculated using Eq. (39), at the convergence of the models’ calibration procedure. All models lead to a very 


Table 5 
Damage models parameters identified using DP600 steel. 

Model Damage Parameters 
Parameters DP600 steel 

B-W Dy 0.9653, 
Dz 0.0002 
D3 0.9036 
Ds 0.9677 
Ds 1.4999 
De 0.8784 

M-W Dy 0.84 
Dz 2.086 
D3 —3.703 
Ds 0.3577 

X-W Cy 1.5791 
Cy 1.0 
C3 0.9671 
C4 0.8866 

L-HC=1/3 Cy 4.0985 
C2 1.2510 
c3 0.9821 

MMC Cy 0.1 
Co(MPa) 505.956 
c3 0.9639 

H-C a 1.0961 
b 0.9753 
c 0.1053 

Q-K Cy 1.9052 
Cy 1.5006 
C3 3.3479 
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similar value of the average calibration error, lower than 12%, except the M-W. 

Fig. 18 shows the relative calibration error, calculated for each of the 6 tests considered in the calibration (tests #1 to #6 for DP600 
steel, in Table 2), which corresponds to the gap between the experimental fracture strain and the one predicted by the model. For this 
material, all models are unable to capture the ductile fracture for the pure shear test (test #2), combined shear-tension test (test #3), 
and Nakazima test (test #6). The 2D fracture locus obtained with the seven uncoupled ductile fracture models are presented in Fig. 19. 
All models, except M-W, describe the decrease of the fracture strain with the decrease of the triaxiality from —0.6 to —1/3. As pre- 
viously mentioned, this exception noticed for the M-W model is related with the fact that this model imposes a cut-off value for 1 = 
—1/3.All describe the increase of the fracture strain with the decrease of the triaxiality from —1/3 to 0. For the low range of positive 
stress triaxiality (0<7<1/3), the material ductility predicted sharply decreases as a function of the increase of 7 and then it increases 
rapidly until the 7 reaches 1/3. In the high range of positive stress triaxiality (7>1/3), the fracture strain starts to decrease and then, 
increases until the 7 reaches approximately 2/3. In fact, all models except M-W, display an asymptotic behaviour close toy = + 2/3, 
under PSC. Nevertheless, this is not enough to capture the Nakazima test (test #6; 7 = 0.66) fracture strain, leading to errors higher 
than 15% for all models, as shown in Fig. 18. Fig. 19 also highlights that L-H, MMC, H-C and Q-K lead to similar 2D fracture locus for 
n2—1/3 (see also the average calibration error in Fig. 17). 


5.2.2. Analysis of the prediction ability 

For DP600 steel, the experimental tests considered for evaluating the prediction accuracy are the ones from #7 to #18, in Table 2. 
These 12 tests were carried out in the high 7 range. Fig. 20 shows the relative prediction errors of the seven ductile fracture models. It is 
clear that all models are not able to predict the fracture strain for the tensile test performed on a Grooved flat, with R = 2 mm (test 
#13), and with R = 4 mm (test #12), and for the smooth flat subjected to tension (test #7). Moreover, all models lead to similar relative 
error values for all the tests with € = 0. On the other hand, the X-W, L-H and Q-K models present also a similar level of prediction 
accuracy for the tests with € = + 1. 

The 3D fracture loci obtained for the DP600 steel are presented in Fig. 21. In this figure, the 2D fracture locus is displayed by the 
thick dash line, whereas the normalized uniaxial compression, plane strain, and uniaxial tension conditions are represented by the red, 
blue and black-coloured solid lines, respectively. The six experimental tests used for the calibration procedure, from #1 to #6, are 
marked with white-coloured circles. The tests used for the prediction accuracy, from #7 to #18, are displayed with the red-coloured 
circle marks. The fracture strain predicted for each test condition (y, €) defined in Table 2 is marked with blue-coloured circles. The 
relative errors for each test (#1 to #18) are denoted by black thin lines. For this material, the B-W model leads to a quite flat surface 
while the ones obtained by MMC and H-C are quite similar. Nevertheless, for the high range of 7 the surfaces present a quite similar 
trend for all models. Therefore, Fig. 21 shows that all ductile fracture models, except the B-W, underestimate the fracture strain of the 
experimental tests considered in the evaluation of the prediction ability. 


5.3. Comparative analysis of the models’ calibration and prediction tendencies for the two materials 


As previously mentioned, the set of experimental tests considered for the DP600 presents a range for the 7 identical to the one of 
Group 1 for the 2024-T351, but a wider range for €. On the other hand, when compared with Group 2 for the 2024-T351, the range for é 
is identical while the range for 7 is narrower. 

Comparing the average calibration errors obtained for 2024-T351 with Group 1 (Fig. 8, grey bars) and DP600 (Fig. 17), it is shown 
that they are identical for both materials (lower than 14%), meaning that the fracture strain values are globally well described for both 
materials. Nevertheless, for the 2024-T351, a smaller average calibration error is obtained, particularity for B-W, X-W and Q-K models. 

If the comparison is made for 2024-T351 calibrated with Group 2 (Fig. 8, yellow bars) the difference between both average 
calibration errors becomes slightly bigger, because it has smaller values than the one obtained with Group 1, except for M-W and Q-K. 
These results indicate that the average calibration error can be improved by considering a wider range foré, but this depends also on the 
type of tests selected, as previously discussed when analysing the influence of the calibration sets. As previously mentioned, for the 
DP600 steel, none of the models are able to accurately capture the fracture strain for the pure shear test (test #2) and the combined 
shear-tension test (test #3), which present a conflicting trend. The same effect is observed for the 2024-T351 calibrated with Group 3, 
when compared with Group 2. 

Regarding the 2D fracture locus obtained for 2024-T351 (Fig. 8 (a) and (b)) and DP600 steel (Fig. 19), the results show that the 
predicted equivalent strains at fracture for the 2024-T351 aluminium alloy are much lower than those obtained for DP600 steel, as 
expected from the experimental results shown in Table 1 and Table 2. It is also noticed that the fracture models present quite different 
trends for both materials, particularity in the low range of positive stress triaxiality (0<7<1/3), when the comparison is made with 
Group 1, for the 2024-T351. On the other hand, the trend becomes more similar when considering Group 2, for the 2024-T351 (Fig. 8 
(b)), highlighting the impact of the fracture strain considered for close to 1/3. In contrast to the fracture locus obtained for 2024-T351 
Group 1 (Fig. 8 (a)), for DP600 steel (Fig. 19), the fracture strain predicted by all the models exhibits an inflection point for a 7 close to 
0.2. This behaviour was also observed for the 2024-T351, for some models, when using Group 3 (see Fig. 8 (c)), which once again 
highlights the importance of the experimental tests selected. Moreover, the X-W and B-W models also present quite different trends for 
n = —1/3 (see Fig. 8 and Fig. 19), for both materials. For DP600 steel, the fracture strains obtained with X-W and B-W models are lower 
than the ones predicted by the other models, for—1/3<n<0. However, for 2024-T351 calibrated using Group 1, the fracture locus 
predicted by these two models are much higher than the ones corresponding to L-H, H-C, and Q-K models, particularly in the negative 
range of 7. The differences observed between the fracture strains predicted by X-W and B-W models, when compared with the other 
models, can be explained by the lack of experimental data points in this range. This could have an impact on the accuracy of the 
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Fig. 17. Average errors predicted in the calibration procedure for the DP600 steel. 
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Fig. 18. Relative errors of fracture models obtained in the calibration procedure under PSC for DP600 steel. 


calibration procedure, particularity for X-W and B-W models, which are more flexible than the others. Moreover, for 2024-T351 
calibrated using Group 1, the X-W model shows high magnitude of the fracture strain for 7 = +1/3 (zone B in Fig. 8 (a)). This is 
likely a result from the elliptic function used to describe the drop of material ductility due to the influence of £, but also a consequence 
of the high value of fracture strain associated to test #14. 

The percentage average errors associated to all the tests used in the prediction accuracy (tests #7 to #18, for DP600 steel and #18 
to #29, for 2024-T351), and the ones obtained from both the calibration and the evaluation of the prediction ability (tests #1 to #18, 
for DP600 steel, and #7, #8, #11, #14, #16, #17 and from #18 to #29, for 2024-T351 calibrated using Group 1; and #2, #6, #8, #10, 
#13, #16 and from #18 to #29, for 2024-T351 calibrated using Group 2) are shown in Fig. 22. From these results, it is interesting to 
note that the X-W model is the only one presenting similar average prediction errors for both materials, when Group 1 is adopted in the 
calibration procedure, which are among the lowest ones. However, the 2D fracture locus are quite different for both materials (see 
Fig. 8 (a), (b) and 19 (c), respectively). Fig. 22 highlights that for the DP600 steel, the B-W and X-W models present more accurate 
prediction than the others models, for the tests from #7 to #18. On the other hand, there are slightly differences in the average error 
values of the seven selected models when considering all tests (see white-coloured bars). For the 2024-T351 aluminium alloy, when 
considering Group 1, the M-W and X-W models present clearly lower values of average prediction errors than the other models (see 
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Fig. 19. Comparison of 2D fracture locus obtained with the seven uncoupled models for DP600 steel. 


DP600 steel 
Relative errors of 12 tests 


oe of oo of 
| oe © © geom i 


300 


) 
i) 
s 


150 


Relative errors (% 


100 


B-W M-W X-W L-H MMC H-C Q-K 


Uncoupled models 


Fig. 20. Relative errors of fracture models prediction for the 12 tests from #7 to #18 for DP600 steel. 


medium grey dashed bars, in Fig. 22). The highest average prediction error values obtained for this material can be explained by the 
fact that all the models, except M-W and X-W, are not able to properly predict the fracture strain for the upper bound assumption (i.e., 
€ = +1), particularly, for smooth and small notched bars subjected to tension (tests #20 and #29) (see Fig. 9 (a)). For DP600 steel, 
Fig. 20 shows that all the models describe slightly better the fracture prediction at the tests corresponding to = +1 than for the 2024- 
T351 aluminium alloy, calibrated with Group 1. Nevertheless, when the comparison is made with the calibration performed using 
Group 2 for the 2024-T351 aluminium alloy, one can see an opposite trend. 

Regarding the 3D fracture locus obtained for 2024-T351 Group 1 (Figs. 10-16), the results show that all models overestimate the 
experimental tests for € = + 1, except for the M-W model. The adopting of Group 2 only changes this trend for the B-W model. 
However, for the DP600 steel, all models, except B-W, underestimate the tests corresponding to € = + 1, as shown in Fig. 21. 


6. Conclusion 


A comparative study on the effect of the type and number of tests used for the calibration of seven ductile fracture models on their 
capability to forecast the failure behaviour of the 2024-T351 aluminium alloy and DP600 steel was presented. The calculated error 
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Fig. 21. 3D fracture loci obtained for the: (a) B-W, (b) M-W, (c) X-W, (d) L-H, (e) MMC, (f) H-C and (g) Q-K models calibrated for the DP600 steel. 


between the predicted and the experimental fracture strains allows to assess both the prediction performance of the calibrated models 
and, consequently, the adequacy of the chosen experimental dataset used for models’ calibration. It is noteworthy, that the available 
experimental data was split in two subsets, for both materials. The first one includes tests that were carried out under PSC. The second 
subset was not used in the models’ calibration, but in the assessment of the ability of the calibrated models to forecast ductile failure. 
For the 2024-T351, the set of tests performed under PSC was divided into three subsets of 6 tests, which covering some chosen stress 
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Fig. 22. Average errors of fracture models for all the tests used in the calibration and prediction for DP600 steel and 2024-T351 aluminium alloy 
(Group 1 and Group 2). 


states. Based on the analysis of the obtained results, the main conclusions are drawn below. 

For 2024-T351 aluminium alloy, it was shown that different types and numbers of tests, lead to distinct calibration results; but low 
average prediction errors are found for the more flexible models. In fact, these models present similar average prediction errors when 
calibrated with 6 or 17 tests, covering a similar range for the variables that characterize the stress state and fracture strain values. It 
was also noticed that low models’ predictivity nearby the specific stress triaxiality values of 7 = —1/3 and n = +1/3 (“angular points” 
or cut-off for M-W model for 7 = —1/3), for which the fracture strain predicted is quite sensitive to the type of tests selected. Therefore, 
the design of specific experimental tests for which the stress triaxiality are nearby these critical values is desirable. 

The calibration of the fracture models using the DP600 steel experimental dataset, which is more restrictive than the one available 
for 2024-T351, was performed in an attempt to generalize the tendencies of the calibration and verification procedures achieved for 
the aluminium alloy. Indeed, it was shown that accurate results in the calibration procedure and in the prediction ability can be 
obtained when the experimental set used in the calibration covers a range as wide as possible of both the stress triaxiality and the Lode 
parameter. The range covered for the Lode parameter is particularly critical for the more flexible models. For the materials under 
analysis, the X-W model seems to be the one leading to a more accurate description of all tests. 

It was shown that an experimental set of at least 6 tests is sufficient for an accurate calibration of the most flexible uncoupled 
models and similar calibration and verification results may be achieved, with what it is obtained using 17 tests. This involves choosing 
experimental tests for the calibration procedure that cover a wide range of stress states. Thence, when the calibration procedure is 
applied using a similar range for 7 and a wider range for €, the average calibration errors obtained for all models tend to reduce. Last 
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but not least, the use of some experimental tests characterized by identical values of stress triaxiality and Lode parameter, but different 
fracture strains, can be detrimental for the calibration procedure. The origin for this discrepancy in the experimental fracture strain 
value can be related, among other sources, with the uncertainty on the location and measurement of the principal strains, the quality of 
the machined samples, the heterogeneity and scattering of the material properties, etc. In this context, the use of the same type of 
specimens (e.g. plate or tubular) can help reducing some sources of scattering. Nevertheless, it will be profitable to use probabilistic 
approaches that can handle different sources of scattering to develop robust calibration procedures. 
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